Quantitative analysis of the dripping and jetting regimes in co-flowing capillary jets 



O 

(N 

(N 

>> 
i 

q 

O 

•i-H 

C/5 
>Y 

^ : 
& 

(N : 
> 

oo : 

(N . 
(N ■ 



Maria Luisa Cordero 1 , Frangois Gallaire 2 and Charles N. Baroud 1 
1 LadHyX and Department of Mechanics, Ecole Polytechnique, CNRS, 91128 Palaiseau, France 
Laboratory of Fluid Mechanics and Instabilities (LFMI), 
Ecole Polytechnique Federale de Lausanne (EPFL) EPFL /STI/IGM/LFMI, CH-1015 Lausanne, Switzerland 

We study a liquid jet that breaks up into drops in an external co-flowing liquid inside a con- 
fining microfluidic geometry. The jet breakup can occur right after the nozzle in a phenomenon 
named dripping or through the generation of a liquid jet that breaks up a long distance from the 
nozzle, which is called jetting. Traditionally, these two regimes have been considered to reflect the 
existence of two kinds of spatiotemporal instabilities of a fluid jet, the dripping regime correspond- 
ing to an absolutely unstable jet and the jetting regime to a convectively unstable jet. Here, we 
present quantitative measurements of the dripping and jetting regimes, both in an unforced and 
a forced state, and compare these measurements with recent theoretical studies of spatiotemporal 
instability of a confined liquid jet in a co-flowing liquid. In the unforced state, the frequency of 
oscillation and breakup of the liquid jet is measured and compared to the theoretical predictions. 
The dominant frequency of the jet oscillations as a function of the inner flow rate agrees qualitatively 
with the theoretical predictions in the jetting regime but not in the dripping regime. In the forced 
state, achieved with periodic laser heating, the dripping regime is found to be insensitive to the 
perturbation and the frequency of drop formation remains unaltered. The jetting regime, on the 
contrary, amplifies the externally imposed frequency, which translates into the formation of drops 
at the frequency imposed by the external forcing. In conclusion, the dripping and jetting regimes 
are found to exhibit the main features of absolutely and convectively unstable flows respectively, 
but the frequency selection in the dripping regime is not ruled by the absolute frequency predicted 
by the stability analysis. 



I. INTRODUCTION 

A liquid column that exits from a nozzle in the form of 
a jet tends to break up into droplets due to surface ten- 
sion. During the breakup process other forces are also 
present and oppose the action of surface tension. In par- 
ticular, viscous forces damp the growth of the varicose 
deformations of the jet that lead to the pinch off and iner- 
tial forces promote the formation of a long liquid thread. 
Gravitational effects can also oppose the action of sur- 
face tension by elongating a falling liquid thread. These 
forces define the parametric space which determines the 
breakup dynamics of the liquid jet. 

Liquid jets exhibit a large range of breaking patterns in 
different regions of the parameter space [l[ . As a general 
rule, the system exhibits a regime of periodic drop emis- 
sion near the nozzle when capillary forces are dominant 
over other forces. This phenomenon is named dripping 
[Fig. [Ha)] and produces monodisperse drops whose size 
is typically determined by the nozzle size. As inertial or 
viscous forces become comparable to capillary forces the 
system transitions to a regime characterized by the for- 
mation of a long jet. Its radius undulates and it breaks up 
into drops far downstream from the nozzle. This regime 
is known as jetting [Fig. QJ b)] and usually exhibits a high 
dispersion in the drop size. 

The transition from dripping to jetting in a jet that 
falls in a quiescent outer fluid under the effect of grav- 
ity has been extensively studied, for example as a model 
of chaotic dynamics On the scale of microfluidics, 

early studies of the transition were performed by Gahan- 
Calvo 0]. In comparison, the dripping to jetting tran- 



sition in confined liquid jets that flow in an immiscible 
carrier fluid has received less attention but interest in this 
geometry has increased since the recent development of 
microfluidic devices for the generation of emulsions Q. 
One of these devices is the so-called co-flowing geome- 
try [9J, in which the inner fluid is injected into a carrier 
fluid that is flowing in parallel [Fig^a)]. Gravitational 
effects are typically negligible in these systems. 

Recently, Utada et al. .9] have studied this system 
when the transversal dimensions of the container chan- 
nel is large compared to the jet radius and they have ob- 
served a dripping to jetting transition when viscous or in- 
ertial forces become comparable to capillary forces. This 
study was followed by several that considered the case 
of a highly confined jet with negligible inertia, which we 
concentrate on below. As shown by Guillot et al. [Tol [TTI| . 
the transition in this case is ruled by the capillary num- 
ber Ca, which measures the ratio of viscous to capillary 
forces and whose critical value depends on the value of 
the confinement. 

The theoretical efforts to understand the dripping to 
jetting transition in these devices have mainly concen- 
trated on the linear stability of the liquid jet. Indeed, a 
liquid jet is unstable to axisymmetric perturbations of its 
base columnar state [l2[. In the case of a flowing jet, per- 
turbations can be modeled as exponential functions with 
complex frequencies and wavelengths. These perturba- 
tions can grow both in time and space while traveling 
along the jet interface [Fig. HJb)]. If the growing per- 
turbations travel both downstream and upstream, the 
instability is said to be absolute. In analogy with the 
high Reynolds number jets d, [IH, [ijj , this case has been 
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Figure 1: Snapshots of the system in the dripping (a) and in the jetting (b) regimes. Experimental conditions are (a) Qi = 
5 /xL/min, Q e = 30 /xL/min and (b) Qi = 14 /iL/min, Q e = 30 /xL/min. 
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Figure 2: (a) Co-flowing geometry. Inner fluid has viscosity rji 
and is injected at flow rate Qi, while outer fluid has viscosity 
r/ e and is pushed at flow rate Q e . (b) Schematic picture of the 
propagation of perturbations in a jet whose interface moves 
with velocity U. v~ and v + represent the edge velocity of a 
wave packet composed of unstable perturbations. 



related to the dripping regime, since instabilities invade 
the whole jet and particularly the exit of the injection 
nozzle. Conversely, if all growing perturbations travel 
downstream the instability is said to be convective and 
this has been related to the jetting regime, in which the 
fluid velocity is high enough to advect the instabilities 
and induce breakup away from the nozzle. 

Guillot et al. 0, E| and Herrada et al. [HI have per- 
formed stability analyses for the confined liquid jet in a 
co-flowing liquid. They predict a transition from an ab- 
solute to a convective instability at a critical value of Ca, 
in the limit of negligible inertial forces. These authors 
find a good agreement between this theoretical critical 
value with the observed value of Ca at which the sys- 
tem switches from dripping to jetting, confirming that 
the two transitions are related. However, there exist no 
detailed comparisons between the theoretical predictions 
and the experiments beyond the threshold of the transi- 
tion, which is the subject of this paper. In particular, we 
perform a frequency and wavelength analysis in both the 
jetting and dripping regimes and compare the measure- 
ments with the theoretical predictions. We also demon- 
strate the possibility to force the breakup of the liquid 
jet with an external stimulus, which allows us to control 
the drop size in the jetting regime. 



In the next section we present the setup used for the 
experiments. Section Mil presents the experimental re- 
sults. In Section IIII Al we present the measurements of 
the flow in the absence of external forcing. Then, in Sec- 
tion IIIIBl we study the response of the flow to an external 
perturbation induced with laser heating. In Section IIVI 
our experimental results are compared with the predic- 
tions of existent theories concerning confined liquid jets. 
Finally, Section [V] summarizes the results and conclu- 
sions. 



II. METHODS 

Our experiments are performed inside a straight 5 cm 
long microchannel of square cross section. The mi- 
crochannel mold is produced by lithography using four 
layers of a photoresist of thickness 50 fj,m each (Etcrtec 
XP-800-20) on a glass slide. The channel is then repli- 
cated in PDMS prepared in a 10:1.5 ratio of PDMS base 
with curing agent. Three inlets and one outlet are drilled 
in the PDMS block (see Fig. [3]). A glass cylindrical cap- 
illary of external diameter 200 /im and internal diameter 
100 /jm is placed between the first and the third inlets to 
form the co-flowing geometry where the drops and jets 
are produced. The alignment of the glass capillary inside 
the square channel is ensured by the match between the 
size of the channel and the external diameter of the glass 
capillary. After bonding the channel to a glass slide, a 
drop of PDMS, prepared in a 10:0.5 ratio of PDMS base 
with curing agent, is allowed to enter the channel through 
the second inlet. After flowing by capillarity through the 
space between the PDMS walls and the glass capillary 
in both directions, the PDMS reaches a wider section of 
the channel where it stops flowing due to the change of 
Laplace pressure [ll| . The channel is then baked to cure 
the PDMS and ensure the sealing between the glass cap- 
illary to the external channel. In this way, the inner fluid 
is forced to flow inside the glass capillary. 

Fluids are injected at constant flow rates using sy- 
ringe pumps (Kd Scientific 101, Harvard Apparatus 
PHD2000W). The inner fluid is a 80% glycerol and 20% 
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water mixture, whose viscosity and density at 20°C are 
rji = 60.1 mPa s and pi — 1208.5 kg/m 3 respectively. The 
outer fluid is silicone oil (Rhodorsil 47 V 100) of viscosity 
i] e = 100 mPa s and density p e — 965 kg/m 3 . Surface 
tension between both fluids is measured as a function of 
temperature at j(T) = 47 — 0.3 T mN/m, where T is the 
temperature, in degrees Celsius. 

The channel is mounted on an inverted microscope 
(Nikon ECLIPSE TE2000-U) and observed through the 
microscope objective (Nikon, 4x/0.1 NA). Image se- 
quences are recorded at 1000 frames per second with 
a fast camera (Photron Fastcaml024 PCI). An infrared 
laser of wavelength 1480 nm (Fitel Furukawa FOL1425) 
is focused inside the channel through the same objec- 
tive, a distance 215 /im from the nozzle. The power 
of the laser beam is modulated with a signal generator 
(Wavetek model 80) to a sinusoidal signal of amplitude 
P up to 55 mW and frequency /i aser ranging between 10 
and 150 Hz. 

The laser radiation is used to heat the inner fluid by 
absorption. This produces a periodic forcing of the flow 
by locally changing the capillary number Ca = rjiU m /j 
(where U m is the mean velocity of the flow), which con- 
trols the stability of the jet. Indeed, the temperature 
increase reduces both the interfacial tension between the 
fluids 7 and the inner viscosity rji fl7| . Figure 0] shows 
the effect of the temperature increase above room tem- 
perature (25°C) on Ca for a typical value of the mean 
flow velocity U m = 20 mm/s. We find that Ca decreases 
when the temperature increases, which means the laser 
acts by locally destabilizing the flow. 

The temperature increase depends on the laser power. 
For the powers used here we expect the heating to range 
between 5°C and 30°C, based on our previous results |18| . 
which corresponds to a local decrease of Ca between 
25% and 70% approximately. Compared to our previous 
work [l9j], the laser heating here is weaker and the flow 
velocity and fluid viscosity larger, so the laser produces 
a small perturbation that cannot block the flow. 



Figure 4: Capillary number as a function of temperature in- 
crease above room temperature for U m = 20 mm/s. 



III. EXPERIMENTAL RESULTS 
A. Unforced state 

First, we study the flow behavior in the absence of the 
laser heating. The external flow rate is fixed at Q e — 
30 /xL/min, while the internal flow rate Qi is increased 
from 5 /xL/min to 32 /iL/min with a step of 1 /iL/min. A 
transition from dripping [Fig. QJa)] to jetting [Fig. [TJb)] 
is observed at Qi = 14uL/min, similar to what was 
described elsewhere [ll|, HO] • 

The length of the inner fluid column, , is measured 
as a function of Qi and is shown in Fig. [5] The error bars 
are calculated by considering the minimum and maxi- 
mum values measured for during the drop formation 
process. For low Qi, is small and grows only slightly 
with Qi. This corresponds to the dripping regime. At 
Qi = 14 /iL/min, the jet starts to grow, showing that the 
jetting regime has been attained. 
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Figure 5: Length of the inner fluid column as a function of 
the inner flow rate for an external flow rate fixed at Q e = 
30 /xL/min. 

The low dispersion in in the dripping regime shows 
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that drop breakup occurs at a steady position, suggesting 
that it is a regular, periodic process. This contrasts with 
the high dispersion in the jetting regime, which reflects 
the fact that the position of drop pinch off is not as well 
defined and therefore the drop formation is not regular 
in time. 

To confirm this, we analyze the frequency of drop pro- 
duction in both regimes by measuring the frequency of 
the oscillations of the inner fluid column. The thickness 
of the jet is measured as a function of space and time 
from image analysis. This yields the radius of the jet 
r(z, t), where z is the axial direction shown in Fig. HJb) 
and t is time, for each value of Qi. A discrete tempo- 
ral Fourier transform yields the frequency spectrum of 
the oscillations of the jet's radius at each point z of the 
jet. These spectra are averaged in space to obtain one 
spectrum for each value of the inner flow rate. 

The frequency spectra as a function of Qi are shown 
in Fig. [SJa) in gray scale. In the dripping region (Qi < 
14 /iL/min) a sharp peak is observed, whose position in- 
creases linearly with Qi. A typical spectrum for this 
region is shown in Fig. [BJb) . In the jetting regime, 
instead, two regions can be distinguished. First, for 
14 < Qi < 24 /iL/min, a broadening of the peak can 
be observed [see Fig. [6^b)] with a global decrease in its 
position. We call this region "first jetting regime" or Jl. 
Then, for Qi > 24 /iL/min, a sharp peak whose position 
increases linearly with Qi can be observed again. This 
region is called "second jetting regime" or J2. A close 
examination shows that the peak of J2 is associated with 
fluctuations induced by the injection of the inner fluid by 
the syringe pump (see Section IIV|) . 

In contrast to the dripping regime, the jet in the jetting 
regime undergoes undulations characterized by a spatial 
wavelength visible on Fig. [ljb) . The wavelength of these 
oscillations, Ao, could not be accurately obtained from 
discrete spatial Fourier transform due to the insufficient 
spatial resolution of the optical setup. Instead, it was 
deduced from the drop volume Vd and the mean jet ra- 
dius calculated as the temporal and spatial average 
of r(z, f). Assuming volume conservation, we write 

V d = TrR2\ , (1) 

which allows us to calculate Ao- Figure [f^c) shows Ao as 
a function of the inner flow rate. The error bars are cal- 
culated from the dispersion of the drop size distribution. 
In the first jetting regime, Aq increases with Qi and dis- 
persion is large, reflecting the polydispersity of the drops. 
In the second jetting regime, Ao remains constant, with 
little dispersion. 

Finally, the phase velocity of the jet oscillations is cal- 
culated from the frequency and wavelength of oscillations 
as v = Ao/g Xp , where /q XP is the position of the maxi- 
mum of the frequency spectrum for each Qi . As seen on 
Fig. E)Jd), vq remains constant at around 45 mm/s with 
large dispersion for the first jetting regime, due to the 



large error bars of Ao- In the second jetting regime, vq 
increases linearly, with low dispersion. 



B. Forced state 

Next, we study the response of the flow to an exter- 
nal periodic perturbation driven by the laser heating. 
For that, Qi and Q e are kept constant while the laser 
is modulated at different frequencies and powers. With 
the flow in the dripping regime (Q e = 30 /iL/min and 
Qi = 7 /iL/min), no noticeable effect is detected, even 
at the highest laser power. Frequency spectra of the 
jet oscillations are calculated for each value of the laser 
frequency /i aser , as explained above, and are shown in 
Fig.[7ja) for a laser power of 55 mW. These spectra show 
the dominance of the frequency peak at /q XP which cor- 
responds to the unforced case. Only at high laser power 
does the peak corresponding to /j aser appear but its am- 
plitude is much smaller than the dominant one, thus not 
changing noticeably the drop generation. 

The same experiment is performed with the flow in the 
first jetting regime (Q e = 30 /iL/min, Qi — 20 /iL/min). 
For low laser power, no response of the jet is observed. 
When the laser power is increased, however, the jet's os- 
cillations synchronize to the laser frequency for some fre- 
quency ranges. This can be observed in the frequency 
spectra of the jet oscillations, shown in Fig. [7^b) as a 
function of the laser frequency for P = 55 mW (simi- 
lar results were obtained for all laser powers higher than 
11 mW). A strong dominance of the peak correspond- 
ing to the laser frequency can be observed in the range 
/laser e [55, 120] Hz. Note that harmonics of /i aser are 
also present in the spectrum of Fig. [7jb) . In the range 
/laser <- 55 Hz the dominant peak corresponds to the 
first harmonic 2/i aser , and therefore in this region the 
jet's oscillations are also controlled by the laser induced 
perturbations. 

Two main effects are observed in the frequency ranges 
where synchronization occurs: the length of the jet de- 
creases as the laser power increases (not shown) and the 
drop size changes as a function of /i aser - Figure E^a) 
shows snapshots of the jet forced by the laser with power 
P = 55 mW at various frequencies. One can observe 
that the jet shortens when the laser is turned on and 
that the drops are larger than without the laser forcing, 
decreasing with increasing /i aser - 

By analysis of the drop size and considering the mean 
radius of the unperturbed jet, the wavelength of the 
forced jet oscillations A/ is computed using Eq. (JXJ) . Fig- 
ure IHb) shows A/ as a function of the laser frequency. 
In the range /i aser € [55, 120] Hz, the wavelength of 
oscillations decreases with the laser frequency, in close 
agreement with the relation /i aser A/ = vq. Here, vq is 
the phase velocity calculated from the dominant wave- 
length and frequency of the unforced state at the same 
flow rates. Note that for fi^j- < 55 Hz, although the 
data are more scattered, the measured wavelengths are 
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Figure 6: (Color online) (a) Frequency spectrum (in arbitrary units) of the jet oscillations as a function of Qi. The red thin and 
blue thick arrows mark the flow rates Qi — 7 ^L/min and Qi = 20 /^L/min, which correspond to the experimental conditions 
of Figs. EI a) and (b) respectively. The frequency spectra for these two particular cases are shown in (b) and represent a 
typical frequency spectrum in the dripping regime (red thin line) and in the jetting regime (blue thick line), (c) and (d): 
Wavelength and phase velocity of the jet oscillations respectively as a function of Qi. In (a), (c) and (d), dashed vertical lines 
are used to separate the different behaviors: Dripping (D), first jetting regime (Jl) and second jetting regime (J2). In all cases 
Q e = 30 /iL/min. 



consistent with the relation 2/j aser A/ = vq. This again 
confirms the dominance of the first harmonic of /i aser in 
the frequency spectrum. 

Finally the dispersion in the drop size, calculated as 
the standard deviation of the drop volume divided by the 
average, is shown in Fig. [5] for three values of /i aser - Re- 
markably, the dispersion in the drop size, which is large 
in the jetting regime, decreases by a factor of about 5 in 
the presence of the laser forcing. 



IV. COMPARISON BETWEEN THEORY AND 
EXPERIMENTS 

Herrada et al. [HI and Guillot et al. [HI [HI have ad- 
dressed the stability problem of the confined, infinitely 
long column of liquid shown in Fig. [2][a), when iner- 
tial effects are negligible. While Herrada et al. solve 
the problem analytically, Guillot et al. use, in addition, 
the lubrication approximation. Both theories, therefore, 
present similar results, although the latter is expected 
to be less accurate for perturbations with small wave- 



lengths compared with the jet radius. In the parameter 
space explored in our study, however, this difference is 
irrelevant since the measured wavelengths are at least 10 
times larger than the jet radius [Fig. Etc)]. 

The analyses of Herrada et al. and Guillot et al. yield 
a dispersion relation that relates the complex wave fre- 
quency uj to the complex wave number k in the form: 

u = vk + ia(k), (2) 

where v is proportional to the velocity of the jet interface 
U, which depends on the flow rates and viscosities of the 
fluids. The growth rate cr(fc) is a nonlinear function of 
the wavenumber k, the fluid viscosities, flow rates and of 
the capillary number Ca. As a consequence of the nonlin- 
earity the propagation of perturbations occurs in a dis- 
persive way, meaning that perturbations with different 
wave numbers travel with different velocities along the 
jet surface. In particular, unstable perturbations travel 
with velocities ranging between v~ and v + , schematically 
shown in Fig. &h). The sign of v~ determines if the flow 
is convectively or absolutely unstable. If v~ < 0, then 
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Figure 7: Frequency spectrum of the jet oscillations as a 
function of /i aser for (a) a dripping flow (Q e — 30 /iL/min, 
Qi = 7 /tL/min) and (b) a jetting flow (Q e — 30 /iL/min, 
Qi = 20 /iL/min) forced with laser heating. In both cases 
P = 55 mW. The dashed line in (b) shows the position of the 
dominant frequency f ^ in the unforced case. 
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Figure 8: (a) Images of the jet perturbed with the laser heat- 
ing at a constant power P = 55 mW and different frequencies. 
The scale bar corresponds to 500 /im. The vertical dashed 
line marks the position of the laser focus, (b) Wavelength of 
the jet's radius oscillations as a function of the laser frequency 
•/"laser ^ or different laser powers. The solid line corresponds to 
the dominant wavelength of the unperturbed jet's oscillations. 
The dot-dashed line corresponds to the curve /jaggj-A = vq, 
while the dashed line corresponds to 2/j ager A = wo- Here, 
Qi = 20 /iL/min and Q e = 30 /iL/min. 



there are unstable perturbations that travel upstream, 
making the flow absolutely unstable. On the contrary, if 
v~ > all unstable perturbations travel downstream and 
the flow is convectively unstable [HI, [H| . Herrada et al. 
and Guillot et al. have found the transition between an 
absolutely and a convectively unstable flow (A/C transi- 
tion) to depend on the capillary number Ca and the A/C 
transition in capillary jets has been traditionally related 
to the transition between dripping and jetting [14| . 

In the following, we test the predictive power of the 
theories by extracting several quantities that can be com- 
pared to our measurements. For that we note first that, 
contrary to the hypothesis of Herrada et al. and Guil- 
lot et al., our confining channel has square cross section 
instead of being cylindrical. As shown also by Guil- 
lot et al. this mainly affects the theoretical results 
when the confinement of the jet 2Ri/h becomes close 
to one. We measure the confinement in our experi- 
ments in the range [0.38, 0.5] and therefore the geom- 
etry of our confinement channel can be considered, for 
calculation purposes, as a cylinder of equivalent radius 
R e = hj ^/tt = 113 (Mm. 

When our experimental parameters (A = 0.6, 7 = 



40 mN/m, Q e = 30 /iL/min and R e = 113 /zm) are in- 
jected into the theories considered here, the A/C transi- 
tion is computed at Qi = 26.2 ^L/min with Guillot's the- 
ory and Qi — 23 /iL/min with Herrada's theory. These 
values are of the same order of magnitude than our mea- 
surement at Qi = 14 /iL/min. The disagreement in the 
exact value of the critical external flow rate could be 
related to the geometric differences in cross section be- 
tween our square channel and the circular assumption in 
the theory. 

Next, we compare the frequency behavior observed in 
the unforced flow regimes to the theoretical predictions. 
According to the spatiotemporal stability analysis, the 
behavior of absolutely and convectively unstable flows 
differs strongly. Absolutely unstable flows are dominated 
by a single, well defined frequency and therefore behave 
as self-sustained resonators pulsating at the absolute fre- 
quency /o- This is because the perturbation with zero 
group velocity is unstable and, since it grows in place, 
it ultimately contaminates the whole domain [22l424l |. In 
convectively unstable flows, conversely, instabilities are 
advected downstream and, without external noise, the 
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Figure 9: Dispersion of the drop size distribution as a function 
of the laser frequency for the experimental condition Qi — 
20 /iL/min and Q e — 30 ^iL/min. 



flow would go back rapidly to the base state. Convec- 
tively unstable flows therefore amplify a broad range of 
unstable frequencies and the most amplified one, also 
noted for simplicity fo, corresponds to the one with high- 
est spatial growth rate [2l( . 

With this in mind we use the dispersion relations from 
Eq. ([2]) to compute the frequency expected to dominate 
the flow, for both theories, as a function of the inner flow 
rate for an outer flow rate fixed at Q e = 30 /iL/min. In 
the absolute region the dominant frequency fo was cal- 
culated as the frequency of the perturbation with zero 
group velocity. In the convective region, /q was com- 
puted as the frequency of the mode with largest spatial 
growth rate — ko.i, which is deduced from an inversion of 
the dispersion relation uj — u>(k) to k = k(ui) with to € 
The result is shown in Fig. [TUJ where the lines represent 
the theoretical predictions and the gray scale image re- 
calls the experimental results of Fig. [6ja) . Both theories 
predict fo to exhibit a sharp decrease for small Qi and 
to grow slightly near the transition. In the convective 
region, fo decreases slowly after the transition and then 
grows again slowly for larger Qi. 

In the following sections we proceed to a detailed com- 
parison between the theoretical predictions and our ex- 
perimental observations in the different flow regimes. 



A. Dripping regime 

Our measurements of the flow in the dripping regime 
reveal most of the features of an absolutely unstable 
flow. Without external forcing, the dripping flow ex- 
hibits one well defined natural frequency that dominates 
the frequency spectrum, as shown in Figs. [SJa) and (b). 
Moreover, when externally forced with the laser heat- 
ing, the flow in the dripping regime does not respond to 
the forcing and remains dominated by the same natural 
frequency that dominates the unforced state [Fig. EJa)]. 
However, this natural frequency does not appear to be 
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Figure 10: Frequency of the dominant perturbation predicted 
by the theories of Guillot et al. (solid line) and Herrada et al. 
(dashed line) as a function of the inner flow rate normalized by 
the flow rate at which the A/C transition is predicted to occur. 
In this calculation, Q e = 30 /iL/min, r\i — 60 mPa s, r/ e = 
100 mPa s, R e = 113 (im and 7 = 40 mN/m. The underlying 
gray scale present our data of Fig. [6ja) with Qi normalized 
by the observed transitional flow rate Qi = 14 /iL/min. 



linked to the absolute frequency predicted theoretically. 
Indeed, the theories predict fo to decrease steeply as a 
function of Qi for small Qi (Fig. [TU)) . which would result 
in the formation of larger drops at a smaller rate. On the 
contrary, we measure a linear increase in /g X ^ with Qi, 
which corresponds to the formation of drops of constant 
volume at a higher rate. 



B. First jetting regime 

In the first jetting regime, a broad range of frequencies 
are amplified, as one would expect if the flow is con- 
vectively unstable. Moreover, the Jl regime exhibits a 
decrease of /g X ^ with the inner flow rate, which is co- 
herent with the theoretical predictions in the convective 
regime near the A/C transition. 

The hypothesis of a convectively unstable flow in the 
Jl regime is further confirmed by our experiments with 
the laser forcing. These experiments show the synchro- 
nization of the jet oscillations with the laser frequency, 
demonstrating the ability of the jet to amplify externally 
induced perturbations. Furthermore, a threshold laser 
power is necessary to synchronize the undulations of the 
jet to the laser frequency, but beyond this threshold the 
spectral response of the system is independent of the laser 
power [see Fig. HJb)]. This is coherent with the linear 
picture proposed here where only the frequency of the 
external forcing is relevant to perturb the flow. 

Finally, the amplified perturbations are found to fol- 
low the relation /j aser A/ = vq. At leading order, this is 
coherent with the dispersion relation of Eq. @ . 
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C. Second jetting regime 



SUMMARY AND CONCLUSIONS 



In the second jetting regime a single, well defined peak 
dominates the frequency spectrum of the jet oscillations, 
which, at first, contradicts the hypothesis that the jetting 
regime corresponds to a convective instability. However 
the dominant frequency /q X ^ increases linearly in the 
J 2 regime, which allows us to relate this frequency to a 
perturbation associated with the syringe pump used to 
inject the inner fluid. In fact, the stepping frequency of 
the syringe pump increases linearly with Qi as 



4Q» 



TrD 2 dx ' 



(3) 



where D — 4.61 mm is the syringe diameter and dx = 
0.088 fim the syringe pump step. While / pump is too 
high compared to /q X ^ in the J2 regime, its fourth 
sub-harmonic (/ pU mp/4) matches it precisely, as shown 
in Fig. QTJ Note also that the second sub-harmonic 
(/pump/2) is also present in the frequency spectrum of 
the unforced state. 




10 15 20 25 

Q. (nL/min) 



Figure 11: Frequency of the perturbation induced by the in- 
jection of the inner fluid. The solid line represents the pump- 
ing rate [Eq. Q], while the dashed and dot-dashed lines show 
its second and fourth sub-harmonics. The latter is amplified 
by the flow. 

These results suggest that there is also a broad range of 
unstable frequencies in the J2 regime but the frequency 
imposed by the pumping rate of the inner fluid selects the 
dominant one. As can be seen from Fig. II H the transition 
from the Jl regime with no dominant frequency to the J2 
regime, with extrinsic frequency selection, occurs at the 
flow rate where the frequency imposed by the syringe 
pump enters the range of unstable frequencies. These 
results are coherent with a convective nature of the flow 
instability in the second jetting regime. 



The frequency response of the system can be summa- 
rized from Fig.[5Ja). In the dripping regime there is only 
one, well defined frequency which is selected intrinsically, 
although not directly linked to the absolute frequency of 
Refs. |To|, [HI- In the jetting regime, the zone of unsta- 
ble frequencies is wider and with a less marked peak. 
The response of the Jl regime is typical of unselective 
noise amplification, with the whole range of unstable fre- 
quencies present in the spectrum, while the frequency 
selection in the J2 regime is extrinsic and is due to the 
frequency imposed by the syringe pump. 

In conclusion, we have given experimental evidence to 
confirm that the jetting regime in a confined co-flowing 
system is related to a convective instability of the in- 
ner thread. Conversely, other mechanisms are likely in- 
volved in the dripping regime and determine the rate 
of drop formation. The frequency selection in this case 
could be influenced by a mechanism akin to the shear 
induced breakup described by Umbanhowar et al. [25| . 
This would suggest that the influence of the growth of hy- 
drodynamic instabilities of confined co-axial jets in this 
regime is less important than the local effect of shear at 
the nozzle, for our particular geometry. 

Confirmation of the instability nature of the flow in 
the jetting regime is important to understand the mech- 
anisms that control drop formation when capillary and 
viscous effects are in competition. In particular, spa- 
tiotemporal stability analysis predicts high inherent dis- 
persion in drop sizes in the jetting regime. This is be- 
cause noise sources, such as micrometer sized fabrication 
defects of the channel walls or fluctuations in the pressure 
and flow rates, are unavoidably present. While this kind 
of perturbations are rapidly damped and usually negli- 
gible in microfluidics, they are exponentially amplified 
by a convectively unstable flow. Understanding the flow 
nature provides a way to control this polydispersity by 
inducing the amplification of a preferred single frequency 
with a controlled forcing. Here, we do so with periodic, 
localized laser heating, demonstrating the possibility of 
externally and dynamically tuning the frequency of drop 
formation and the drop size with low dispersion. 
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